Could the compact remnant of SN 1987A be a quark star? 
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ABSTRACT 

The standard model for Type II supernovae explosion, confirmed by the de- 
tection of the neutrinos emitted during the supernova explosion, predicts the 
formation of a compact object, usually assumed to be a neutron star. However, 
the lack of the detection of a neutron star or pulsar formed in the SN 1987A still 
remains an unsolved mystery. In this paper we suggest that the newly formed 
neutron star at the center of SN1987A may undergo a phase transition after the 
neutrino trapping time scale (~10 s). Consequently the compact remnant of SN 
1987A may be a strange quark star, which has a softer equation of state than that 
of neutron star matter. Such a phase transition can induce the stellar collapse 
and result in a large amplitude stellar oscillations. We use a three dimensional 
Newtonian hydrodynamic code to study the time evolution of the temperature 
and density at the neutrinosphere. Extremely intense pulsating neutrino fluxes, 
with submillisecond period and with neutrino energy (> 30 MeV) can be emitted 
because the oscillations of the temperature and density are out of phase almost 
180°. If this is true we predict that the current X-ray emission from the compact 
remnant of SN 1987A will be lower than lO^^erg s~^, and it should be a ther- 
mal bremsstrahlung spectrum for a bare strange star with surface temperature 
of around ~ 10^ K. 



Subject headings: supernova 1987 A: phase transitions: dense matter - quark 
stars 
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Introduction 



Supernova 1987A, the brightest supernova seen in modern times, was a milestone 
i n astronomy and fundamental phys ics. The detection of the neutrinos from SN 1987A 
( iBionta et aLlll987l : iHirata et aLlllQSTI ). confirmed the basic predictions of the physical mech- 
anisms for Type II supernovae. The explosion is triggered by the collapse of massive stars. 
The core collapse of the massive star (with mass greater than IOMq) is due to electron cap- 
ture and photodisintegration, and it is halted when the center density exceeds the nuclear 
matter density. As the stellar core collapses, the gravitational energy is released through the 
emission of neutrinos. A small fraction of these neutrinos is absorbed, and the heating of 
the neutrinos drives the supernova explosion. 

The supernova SN 1987A was observed in every band of the electromagnetic spectrum, 
from radio to gamma rays, and hence it has a well-measured bolometric light curve. By 
using both the bolometric light curve and the spectral evolution of its Ha line, the hydro- 
dynamic and time-dependent atmosphere models are consistent with a presupernova radius 
of 35±5-R0, an ejecta mass of 18 ± 1.5Mp,, an explosion energy of (1.50 ± 0.12) x 10^^ erg. 



and a radioactive ^^Ni mass of O.OTGSM© (jUtrobiru 120071 ) . The inferred energy (~ 3 x 10^^ 
ergs), temperature (~ 10 MeV) and decay time (~ 4 s) of the first neutrino burst were not 
inconsist ent with what woul d be expected from the production of a neutron star through core 
collapse (IHirata et al.lll987l ). Presently, about 35 plausible, or at least possible, associations 
of pulsars with supernova remnants are known, and these are all relatively young pulsars, 
mostly with Tp = P/2P < 10^ years, where P is the first derivative of the pulsar period P 
( iManchesterl 120071 ) . Almost immediately after the observation of S N 1987A, opti c al att empts 
were made to identify the compact remnant of the explosion (see IManchesterl ( 120071 ) for a 
recent review of the optical and X -ray observations on the supernova). The optical searches 
were done bv lPercival et al.l (Il995r). using the High Spee d Photometer System on the Hubble 
Space Telescope, and by iManchester fc PetersonI (119961 ) . using the 3.9 m Anglo-Australian 
Telescope, respectively. No significant pulsations were observed in the period range of 0.2 
ms to 10 s, with an upper limit for the pulsed emission equivalent to a V magnitude of 
about 25. However, most pulsars are detected at radio wavelengths. An extensive search 
was carried out recently at the Parkes 6 4 m radio telescop e in 2006, at frequencies between 
1390 MHz and 8370 MHz, respectively (IManchesterl l2007l ) . No significant candidate with a 
signal to noise ratio greater than 9.0 was observed, and there is no observat ional evidence for 



the presence of a central point (or ne ar point) source at any wavelength (jPark et al.l 12004 



Graves et al.l l2005l : IShtykovskiyI |2005| ) . The limit on the luminosit y of any point sou rce at 
the center of SN 1987A is of the order of 10^^ erg s"^- 10^^ erg s"^ jManchesterl 120071 ) . 



Therefore, the question of why no compact object is observed at the center of the SN 
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1987A explosion is a rather intriguing one. Several explanations have been advanced to 
explain the lack of a ne utron star/ pulsar. The most obvious one is that the pulsar may not 
be beamed towards us (IManchesterl 120071 ) . The beaming fraction (fraction of the celestial 
sphere swept over by the beam as the star rotates) is not very well determined, but in 
radio band a b eaming fraction o f 20% is usually assumed, and for young pulsars it could 
be even larger (IManchesterl 120071 ) . However, even if the pulsar were not directed along our 
line of sight, its radiation would heat the surrounding supernova remnant, and add to its 
bolometric luminosity. At present, all of the luminosity of SN 1987A can be ac counted for 



by a r adioactive dec ay model consistent with the production of O.O75M0 of Ni (iFryer et al. 



19991 : lUtrobinI 120071 ). It would be also possible that the pulsar magnetic field develops in 
a few decades, and in this case a rapidly spinning neutron star would be still undetectable 
(IManchesterl 120071 ). However, most of the current models assume that the magnetic field 
i s either frozen in th e star, or it is generated during the collapse by the dynamo process 
(IBonanno et al.ll2005l ). Another possibility is that the neutron star furth er collapsed into 
a bla ck hole, whose accretion luminosity is below the observational limits (IBethe fc Brown 
19951 ). Although the duration of the neutrino burst of SN 1987A precludes further collapse 
in the first 20 s, the neutron star could have collapsed eventually into a black hole because of 
either accretion of a s ufficient fall-back mass, or due to some changes in the equation of state 
of the dense matter (IFryer et al.lll999l ). The analysis of the collapse of the initially formed 
neutron star in SN 1987A puts severe constraints on the equation of state of the dense 
matter. However, to explain the collapse of a low mass neutron star an exotic equation 
of state of the nuclea r matter is needed, li ke, for example, an equation of state softened 
by pion condensation (IBethe fc Brownlll995l ). which would contradict other observations on 
supernova remnants. Initially, it is likely that neutron stars formed in Type II supernova 
explosions are obscured by the late-time fallback. Much of this fallback is quickly accreted via 
neutrino cooling, but some material may remain on the neutron star, forming an atmosphere 
that slowly accretes through photon emission. If the neutron star has either a low magnetic 
field, or a low rotational sp in frequency, then the neutron star remnant of SN 1987A can not 
be seen (jPryer et al.lll999l ). 



When studying the neutrino data from the SN 1987A explosion a number of unexpecte d 
features have been found ( Hirata et al. 1987 : Bionta et al.lll987 : Costantini et al.ll2004 . 2007). 
For example, the an gular distribution of the two events, one seen at the K amiokande-II (Kll) 
(IHirata et al.lll987l ) and the other at Irvine-Michigan-Brookhaven (1MB) (IBionta et al.lll987l ) 
are more forward directed than expected. The average cosines of the polar angles for these 
two events are (cos^^^^) ~ 0.3, while (cos^''^*^^) ~ 0.5. Moreover, the energy distribution of 
these two detectors are not in perfect agreement. The mean energy of the neutrinos detected 
at KII is around half of the energy of the neutrinos detected at those detected at 1MB, with a 
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mean energy i^E^^'^^^ ^ 30 MeV. The time distribution of the two events is also very different, 
on average the neutr inos detected by 1MB came ~ 5 s after those detected by KII (see 
Alekseev et aLl (119881 ) for a summary of neutrino arrival time for various detectors). When 
fitted with thermal spectra, the two independent detections do not seem to agree with either 
each other or typical theoretical expectations. Using parameter-free inferential statistical 
methods, it can be shown that the combined KII and 1MB data can be best explained by a 
spectral shape that is enhanced both at the peak and t he tail of the spectrum, a nd depressed 
in between, as compared to the Fermi-Dirac spectrum (jYiiksel &: Beaconj|2007l ). Since these 
methods make no a priori assumptions and do not rely on parameter estimation, they allow 
for more efficient processing of small data samples. While the supernova neutrino spectra 
are expected to be quasi-thermal, modifications due to non-standard effects, like neutrino 
mixing among various flavors, neutrino decay, neutrino- neutrino interactions and/or any 
novel mechanism due to unknown physics may produce a time-integrate d spectrum that 
deviates signiflcantly from a quasi-thermal shape (jYiiksel fc Beaconj 120071 ). 



It is the purpose of the present paper to propose some explanations for the lack of 
observational evidence for the presence of a point source at the center of SN1987 A, as well 
as the apparent discrepancy between the two neutrino detections. 



Since IWittenI (119841 ). following early proposals by lltohl (Il970l ) and iBodmerl (Il97ll ). sug- 
gested that strange quark matter, consisting of u-, d- and s-quarks is energetically the most 
favorable state of the matter, the problem of the existence of strange quark stars has been 
intensively investigated in the physical and astrophysical literature. The possibility that 
some compact objects could b e strange stars remains an interesting and intriguing, but still 
open question. IWittenI (Il984j ) also proposed two ways of formation of strange matter: the 
quark-hadron phase transition in the early universe and conversion of neutron stars into 
strange ones at ultrahigh densities. In the theories of strong interaction quark bag models 
suppose that breaking of physical vacuum takes place inside hadrons. As a result vacuum 
energy densities inside and outside a hadron become essentially different, and the vacuum 
pressure on the bag wall equihbrates the pressure of quarks, thus stabilizing the system. 
If the hypothesis of the quark matter is true, then so me of the neutron stars could ac- 



tually be strange stars, buil t entirely of strange matter (IHaensel et al.lll986l : lAlcock et al. 



1986; Alcock fc Ohnto 



198 



However, there are ge neral arguments against the existence 



of strange stars, e. g. ICaldwell fc FriedmannI (1199 
prope rties and the physics of phase transitions see 



For a general rev iew of strange star 



Cheng et al.l (Il998al ) and iHaensel et al. 



A possibility for the formation of the quark stars is that some neutron stars in low- 
mass X-ray binaries can accrete sufficient mass to undergo a phase transition to become 
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strange stars (jCheng fc Dailll996l ). This mechanism has also been proposed as a source of 
radiation emissio n for cosmological 7 - ray bursts, soft gamma-ray repeaters or othe r astro- 
physical objects (ICheng fc Dailll998bl : ICheng et al.lll998l : iPaczynski fc Haenselll2005l ). Some 
basic properties of str ange stars like mass, radius, cooling, collapse and surface radiation 



have been also studied ( Christiansen et al. 2000 : :ICheng fc Harkoll2003l:Earko &: Chenell 



2000 



2002 



Zdunik et al. 2001: iGondek-Rosinska et al. 



2003 



Bejger fc Haensel 



2004 



Zdunik et al. 



20061 ). Quark stars are expected to form during the collapse of the core of a massive star 
after the supernova explosio n, as a result of a first or se cond order phase transition, resultiri g 
in deconfined quark matter ( Takahara fc Sato 1985al lbl: ICentile et al. 1993 : Dai et al. 1995 ). 
The proto- neutron star core, or the neutron star core, is a favor a.ble en vironment for the 
conversion of ordinary matter to strange quark matter (jchen et al. 2007 ). 



This paper is organized as follows. In Section 2, we discuss various possible phase 
transition processes in a neutron star, and to estimate the characteristic time scale for the 
phase transition. We suggest that the phase transition from neutron star matter to strange 
matter in the core of neutron star may be more favorable. In Section 3, we present the 
numerical simulation results of a phase-induced collapse neutron star. In particular, we show 
why high intensity and high energy neutrinos can be emitted after the phase transition. In 
Section 4, we summarize the energy dissipation processes of bare strange stars and their 
cooling. In Section 5, we give a brief discussion about how the phase transition may take 
place in the compact remnant of 1987A and final remarks are presented. 



2. Phase transitions in high density neutron matter 

During the formation and evolution of neutron stars new states of matter may form 
inside the stars as a result of a phase transition, which may be triggered by the accretion of 
matter, pulsar spindown, collapse of the core of a proto-neutron star in supernova explosions 
etc. In the following we shall restrict our discussion only to phase transitions of the first 
order. As a general physical model we will discuss the nucleation of a phase B in the 
metastable phase A. The nucleation is concerned with fiuctuations of parameters, such as 
the local density or the number of particles i n a metastable drop of phase B which triggers 



the phase transition (ILandau fc Lifshitzlll980l ) . There are two extreme cases of the nucleation 
theory. In the classical regime, the temperature is assumed to be sufficiently high to trigger 
the phase transition by thermal fiuctuation. In the quantum regime, below a characteristic 
temperature, thermal fiuctuations are negligible as compared to the quantum ones. Hence 
quantum fiuctuations can initiate a phase transition via the quantum tunneling effect. The 
relevant thermodynamic potential for the description of phase transitions is the Gibbs free 



- 6 - 



energy (IHaensel fc Schaefferlll982l ). In the following we denote by AF the excess free energy 
of a critical droplet. 

The probability of a local formation of a droplet of phase B via a fluctuation decreases 
very strongly with the increase of the number of baryons A^rop in the droplet. For small 
Adrop, however, the positive contribution of the surface energy to AF prevails over the gain 
in the bulk binding, AF (Adrop) > 0, which makes the droplet unstable with respect to the 
reconversion to phase A. However, at some value A^rop = ^crit, the energy excess due to 
the droplet formation vanishes, AF (Acru) = 0. Therefore the droplets with A^rop > Acrit 
grow spontaneousl y, destabilizing the metastable phase A and inducing the phase transition 
harko et aDl2004h . 

By using the thermod ynamical formalism of nucleation ( Landau Sz Lifshit3 1980 : Haensel &: Schaeffer 



19821 : iHaensel et al.l 120071 ) one can study the nucleation of exotic phases (pion condensate, 
kaon condensate and quark matter) in the high density neutron matter. In the following we 
will briefly review each of these processes. 



2.1. Meson condensation 



The nucleation of the pion condensate in the neutron star core and its astrophys i- 
cal implications y y ere st udied by IHaensel fc Schaefferl (119821 ) ; IHaensel fc Proszynskil (119821 ) ; 
Muto &: Tatsumil (Il990l ). respectively, where a theory of metastability of dense neutron mat- 
ter with respect to the first-order phase transition to a pion-condensed state was also an- 
alyzed. The description of the condensation process is based on the idea of the nucleation 
of the pion-condensed phase in the metastable normal matter, through the appearance of 
spontaneously growing droplets of the new phase. Different paths lead ing from the false to 
the tr ue ground state of dense neutron matter have been considered in IHaensel fc Schaeffer 
(Il982h . 



The calculations performed for realistic models of cold neutron matter yield an interval 
of metastability (on the time scale of the age of the universe), which is as large as half of 
that between the baryon density pn, where the metastability starts, and pcrit, where the 
poten tial barrier (without surfa ce effects) between the true and the false ground state van- 
ishes ( IHaensel fc Schaefferlll982l ). In the case of the pion condensation in hot neutron matter 
(T ~ 5 MeV) the region of metastability (on the time scale of the gravitational collapse of 
a massive star) is much narrower than in cold neutron matter. Consequently, additional 
energy release and entropy generation from the first-order phase transition in hot supercom- 
pressed matter are negligible. The value of the A^rit strongly decreases with the growth of 
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the overcompression APover = P — Pn- Thus, the lifetime of an overcompressed sta te de- 
creases rapidly with increasing ^Pover ( iHaensel &: Schaefferlll982l : iHaensel et al.ll2007l ). For 
the case of the pion coridensat ion the surface tension can be approximated by the expression 
derived by iBaym et al.l (Il97ll ). with the nucleus replaced by a pion condensed droplet, and 
the neutron gas replaced by the ordinary neutron star matter. Then, one can obtain the 
condition for the pion condensate in the neutron star core during a time interval equal to the 
present age of the universe, tu = 1-5 x 10^° years. Depending on t he employed model, the 
required overcompr ession is /\PoverlPo = 0.02 or /S.Pover/Po = 0.05 (IMuto fc Tatsumilll990 
Haensel et alll2007h . 



Under typical conditions in the neutron star core (T < 10^ K) the nucleation proceeds 
via the quantum tunneling through the energy barrier. At much higher temperatures, the 
thermal effects increase the nucleation rate through thermally excited droplet states. On the 
other hand, the growth of the temperature increases Pent- For a newly borr i neutron star 
with T > 10^° K, the nucleation proceeds in the classical (thermal) regime ( iHaensel et al. 



20071 ). Due to the softening of the equation of state, pion condensation could lead to a 
significant decrease of the maximum mass and mome nt of inertia allowable for ne utron star 
models, constructed using such an equation of state (IHaensel fc Proszynskilll982l ). 



The formation of a droplet of kaon condensate in a neutron star core is connected with 
the production of strangeness, and nucleation should involve the weak interaction process. 
The transition from npe-type nuclear matter (consisting of neutrons, protons, and electrons) 
to matter containing strangeness, using a W alecka-type mo del, predicting a first-order kaon- 
condensate phase transition, was studied in iNorsenI (120021 ). The free energy of the droplets 
of the kaon-condensed matter, as well as the density, temperature, and the neutrino fraction 
were obtained. The surface tension of the interface between the normal and condensed 



kaon phase was calculated by IChristiansen et al.l (120001 ). using the non-uniform relativistic 
mean field model. In the approximation in which only linear terms in the curvature of the 
surface are kept, the surface contribution to the thermodynamic potential of a spherical 
droplet \s a = as + lod Rdrop-, where as is the surface tension and Oc is the curvature 
coeffic ient. For a small admixt ure of kaon-condensed droplets in a nucleon matter as = 30 
MeV dChristiansen et al.l [20001 ). 



The case of the kaon condensation is drastically different from the case of the deconfine- 
ment transition, since here there is no intermediate zero-strangeness state, which might allow 
for the fast nucleation followed by a slow but smooth growth of the strangeness containing 
fields. Instead, the thermal fiuctuations responsible for nucleation events must directly in- 
volve the weak-interaction processes which produce kaons. The weak interaction processes 
producing strangeness via the reaction e + N ^ Ue + K~ + (where an additional nu- 
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cleon N is needed for momentum conservation) and n —>■ p + K~ are too slow to create a 
critical drople t of kaon condensate from a density fluctuation during the fluctuation lifetime 
( lNorsenll2002l ). Strangeness can be produced at a reasonable rate from thermal kaon-antikaon 
{K~K^) pairs, but this mechanism can operate only at extremely high temperatures typical 
for protoneutron stars. G enerally, a protoneu tron star cools so rapidly that kaon condensate 
has no time to nucleate (IHaensel et al.l 120071 ) . However, as soon as /ig > i^^-, where a;^_ 
is the minimum energy of a s ingle zero momentu m kaon in dense matter, a spontaneous 



formation of kaons is p o ssible (IHaensel et al. 



studied by iMuto et all (119971 ). 



Muto et al 



200"; 



Th e kinetics of kaon co ndensation was 
(|2000aD and iMuto et~aD J2000bh . Kaon poten- 



tials Uk with values only a little below Uk = — 120 MeV would not be compatible with the 
mass of the Hulse- Taylor pulsar, because of a mechanical instability that is initiated by the 
central densities for which the press ure remains consta nt, and the necessary condition for 
stabihty dM/dpc > is not satisfied (lGlendenningll2000l ). The kaon condensation process in 
neutron stars in the framewor k of the Zima n yi-Mo szkowski model in the relativistic mean 
field theory was considered in iDai fc Cheng] (119971 ). Even though hyperons which may in- 
crease the critical condensation density are not included, kaon condensation may not occur 
in stable neutron stars for these classes of mean field theories. The existence of the antikaon 
condensation phase in neutron stars in the frar neworks of the Gle ndenning-Moszkowski and 
Zimanyi-Moszkowski models was reanalyzed in IWang et al.l (120071 ). The results of this anal- 
ysis show that in the very massive (M = 2 ± O.2M0), and high redshift neutron stars, there 
are still some stiff enough equations of state of neutron matter, so that the pure antikaon 
condensation phase and the mixed phase of normal baryons and antikaon condensation can 
still exist. The mass of the neutron star inferred from the neutrino flux of SN 1987A should 
be of the order of 1.5 — I.7M0, which suggest that the kaon condensate may not occur in 
the compact object of the SN1987A. 



2.2. Quark deconfinement 

When quarks seeds are formed in the core of the neutron stars, they will propagate 
through the entire star, and convert it to the new phase. A phase transition occurs between 
the hadronic and quark phase when the pressures and the chemical potentials in the two 
phases are equal, Ph = Pq, fih = fJ^q, where Ph, Ph and Pq, Hq are the pressures and the 
chemical potentials in the hadron and quark phase, respectively. If the transition pressure 
is less than that existing in the supernova core the transition can occur. 

The change from the metastable neutron matter phase to the stable quark phase oc- 
curs as the result of fluctuations in a homogeneous medium, formed of neutrons, in which 
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small quantities of the quark phase (called bubbles or nuclei) are randomly generated. Since 
the process of creation of an interface is energetically unfavorable, it follows that when a 
quark nucleus is below a certain size, it is unstable and disappears again. Surface effects 
disfavor the survival of small bubbles below the radius Rc (called critical size-the nuclei of 
this size are called critical nuclei or bubbles), which is nothing but the value that extrem- 



izes the thermo dy namical work W necessary to create the bubbles (lAlcock fc 01intdll989 



Harko et al.ll2004f). Only nuclei whose size r is above the value Rc are stable, and can survive 
(ILandau fc Lifshitzlll980l ). The nuclei are assumed to be macroscopic objects containing a 
large number of particles (quarks). 

Following the phase transition to the two flavor quark matter, the two flavor quark 
matter will convert into three flavor quark matter through the re actions u + e~ ^ d + v^. , 



u + e ^ s + i/g, and u + d u + respectively (jPai et al.lll995l ). 



Once the quark phase is formed inside the neutron star, it will propagate throughout 
the entire star. The physical mechanisms of the transition from neutron matter to quark 
matter in an astrophysi cal background have been studied within several models. The flrst 
is due to lOlintd (119871 ). who used a non-relativistic diffusion model. As such, this is a 
slow combustion model, with the burning front propagating at a speed of approximately 10 
m/sec. This is determined primarily by the rate at which one of the down quarks inside the 
neutrons is converted, through a weak decay, to a strange quark: d + u — » s + u. The second 



metho d of describing the conversi on process was flrst sugge sted bv iHorvath fc Benyenuto 
(1988), and analyzed in detail by Lugones et al. (1994) and Lugones fc Benvenuto ( 19951 ). 
who modeled the conversion as a detonation. In this case the conversion rate is several orders 
of magnitude faster than that predicted by the slow combustion model. This model is based 
on the relativistic shock waves and combustion theory. But regardless of the way in which 
the transformation occurs, an initial seed of quark matter is needed to start the process. 

The problem of the combustion of neutron matter is closely related to that of shock 
waves. Let us assume that we have an unburnt fluid which converts (by means of a certain 
reaction) in a burnt fluid. The combustion process must be exothermic if it propagates spon- 
taneously to other regions of the fluid. The condition for spontaneous propagation to other 
region of the fluid is E^y^rnt {P, X) < Eunbumt {P, X), where Ebumt {P, X) and Eunbumt {P, X) 
are the energy densities of the respective fluids, both evaluated at the same thermodynamic 
state. In the case of the transition from nuclear to qua rk matter this condition can be re- 
formulated as Enud — 'iPnud > 4-B (ILugoues et al.lll994l ). where En^d and Pnud correspond 
to nuclear matter. If this condition is not fulflUed, the combustion is no longer exother- 
mic, and so it is not possible. At low enough densities, soft EOS verify Enud » Pnud, and 
Enud ~ mnTiBC^ ■ Heucc there exists an absolute lower limit for the combustion to be possible. 
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given by ub ~ AB/ninC^ = 0.25 [-B/(145MeV)^] fm ^ ( Lugones et al.lll994l ). corresponding 
to a transition density pn- of the order of ptr ~ 4.2 x 10""^^ [i?/(145MeV)^] g/cm^. 

However, it is important to point out that the efficient burning of the nuclear matter can 
take place only when the nuclear matter density is sufficiently high, and in order to convert 
the entire neutron star to a quark star a higher density than the minimum one given by is 
required, and there is always a range of densities at which the transition occurs. The density 
for the deconfinement of baryonic matter with a moderately stiff (or stiff) equation of state 
to two-flavor quark matter is near 8pnuc- For a soft equation of state, the deconfinement 
density may be lower. 

Generally, the flow behind the detonation is sonic. Therefore the typical time scale 
for the transition is = R/cs, where R is the radius of the neutron star and Cg is the 
speed of the sound. A simple phenomenological mo del for the evolutio n of the quark phase 
can be obtained by assuming dr/dt = (r — Rc) /rtr (IHarko et al.ll2004j ). which gives for the 
transition time scale Ttr from a microscopic quark nugget to a quark matter distribution of 
a macroscopic size the expression 



Ti,^10-X'^'^6ln-^s 



:i^ 



where Rq is the neutron star radius in units of 10^ cm, Rq ~ 300-Rc (Harko et al. 2004) 
is the initial size of the quark bubble, and Ng is the number of quark seeds inside the core 
of the neutron star, which could be as large as 10"^^ (lida & Sato 1998). However, it was 
found that the conversion process in a h adronic matter at T = always correspond to a 
deflagration, and never to a detonation (jPrago et al.l 120071 ). Hydrodynamical instabilities 
can develop on the front, and a mixed phase of hadrons and quarks could form. Due to 
the formation of the wrinkles, the conversion velocity can signiflcantly be increased, but this 
increase is not sufficient to transform the deflagration into a detonation. In general, one 
could assume that the conversion process can take place in two steps, with a flrst transition 
from hadrons to ungapped (or 2SC) quarks, followed by a second transition in which a CFL 



phase is produced. The time scale for the flrst transition has been estimated by iDrago et al. 



( 120071 ). and it is of the order of 0.1 — Is for the flrst transition in the case of a laminar front, 
and is much more rapid if the hydrodynamical instabilities are taken into account. The 
second transition lasts only 10~^ s, due to the formation of a convective layer. If the two 
processes take place one after the other, it is even possible that the formation of a diquark 
condensate could accelerate the conversion process by developing a convective layer inside 
the hadronic pha se. The neutrino tra pping can play an important role in the hadron - quark 
phase transition (jVidana et al.ll2005l ). The quantum nucleation of a quark matter drop, and 
therefore the conversion of the hadronic matter to quark matter, is strongly inhibited at 
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T = by the presence of neutrinos. However, in the high temperature regime, the dominant 



nucle ation mechanism is thermal nucleation, and not quantum tunnehng (IHaensel et al. 



20071 ). In general, it takes a sub-milliseconds time interval to convert the matter at the core 



of a neutron core into quark matter. 



3. Neutrino emission from a phase-induced collapse neutron star 

The total energy emitted in the form of neutrinos in a supernova explosion can be 
easily estimated from qualitative considerations. The total gravitational energy that can 
be irradiated is Ej, ^ SGM"^ /5R, where M and R are the mass and the radius of the star, 
respectively. By using for the neutron star a mass of the order of M = (1 — 2)Mq, and 
a radius of i? = 20 km (Mq/M)^^^ , one obtains Eh ~ (1 — 5) x 10^^ erg. This amount 
of neutrino energy does not conflict with the observed data. However, it is difficult to 
understand why in comparing with KII, 1MB has detected a neutrino burst with a time 
delay (~ 5s) and with even higher average energy (~ 30MeV). In discussion section, we 
will argue that a phase transition may take place after the core temperature of the newly 
born neutron star is reduced. The characteristic time scale for the cooling of the core is the 
neutrino trapping time, which is of the order of several seconds. In this section, we use a 
three dimensional Newtonian hydrodynamic code to simulate the neutrino emission from a 
phase-induced collapse of neutron star. 



3.1. Description of numerical code 

First we briefly summarize the numerical code used to simulate the collapse of neutron 
star induced by phase-transition. The three-dimensional simulations are based on Newtonian 
hydrodynamics and gravity. The code has been used to study the gravitational wave emission 
from the phase- induced collapse neutron stars (Lin et al. 2006). We refer the reader to Lin 
et al.(2006) for a detailed discussion of the numerical code. 

The system of equations describing the non- viscous Newtonian fluid flow is given by 

|^ + V-(pv) = 0, (2) 

-(p..) + V.(p..v) + - = -p-, (3) 
|L + V-((r + P)v) = -pvV$, (4) 
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where p is the mass density of the fluid, v is the velocity with Cartesian components Vi 
[i = 1, 2, 3), P is the fluid pressure, $ is the Newtonian potential and r is the total energy 
density, r = pe + pv^/2 , and e is the internal energy per unit mass of the fluid. The 
Newtonian potential $ is obtained by solving the Poisson equation, = AnGp. The 

system is completed by specifying an equation of state P = P{p, e). 

The above hydrodynamics equations ©-(jl]) can be rewritten in a flux conservative 
form, which can be solved numerically using quite standard high-resolution shock capturing 
(HRSC) schemes. An HRSC scheme has the ability to resolve discontinuities (e.g. shock 
waves) in the solution. It can also achieve high accuracy in regions where the fluid flow is 
smooth. In our work, we employ the so-called Roe's solver in the simulations (see Lin et al. 
2006). 

It is not known what is the equation of state (EOS) for the neutron star in the remnant 
of supernova 1987A. We can try all possible existing realistic equations of state in our study. 
However, the main purpose of this paper is to demonstrate that the phase-induced collapse 
of a neutron star can emit extremely intense, pulsating and very high energy neutrinos. 
Therefore instead of using realistic equation of state we will use a polytropic equation of 
state for the initial neutron star. In Section 2, we have discussed various possible phase 
transitions for a newly born hot neutron star. Although both pion condensate and kaon 
condensate cannot be ruled out completely, there are no compelling observation evidence 
for the their existence. For simplicity we will focus our study on the phase transition from 
neutron star to strange star in this paper. We use a mixed phase EOS to mimic a strange 
star covered by normal nuclear matter. 

In the following we want to argue that a newly born quark star should be described by 
a mixed phase EOS. In section 2.2, it has been shown that quark seeds can be spontaneously 
generated inside the core of neutron star when the density exceed the critical density (e.g. 
lida & Sato 1998). These quark seeds formed in the neutron star core can mix with the normal 
hadronic matter by Schwarzschild convec tion, and the length scale Ac for the convective 



motion is given by (IWilson &: Maylelll988l ) 



A, 



^" 1'". (5) 



g {—d\nS/dr) 

where 5" is the entropy per baryon in units of Boltzmann's constant, g is the local gravitational 
acceleration in the core, kr = \u/^, where is the neutrino mean free path, rj = X^pi,/3p, 
where pi, is the neutrino energy density and p the matter density. With A,y = 8 cm. 



Piy/p = 0.1, g = 8 X 10^^ cm/s^ and —dlnS/dr = 5 x 10"^ cm~^ (lDai et al.lll995l ). one 
obtains Ac = 3.5 x 10^ cm. The timescale of the mixing can be calculated from t^. = 
271/ y^g {—d In S/dr) = 0.3 ms. Therefore the convection will result in the mixing of the 
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neutron and quark phases in the large region of the star during the phase transition process 



flGlendennindl2000f ). 

The initial equilibrium neutron star before the phase is given by a polytropic EOS 

P = A;o/°, (6) 

where ko and Tq are constants. On the initial time slice, we also need to specify the specific 
internal energy e. For the polytropic EOS, the thermodynamically consistent e is given by 

e = ^P'^-'. (T) 
i — i 

Note that the pressure in Eq. ([6]) can also be written as 

P = (Eo - l)pe. (8) 

We assume that the phase transition take place at t = 0, then we switch the compact 
object from a polytropic EOS to the EOS of a mixed phase quark star, which consists of 
two parts : (i) a mixed phase of quark and nuclear matter in the core at density higher 
than a certain critical value ptr (quark seeds can spontaneously produce everywhere when 
P > Ptr) (ii) a normal nuclear matter region extending from p < ptr to the surface of the 
star. Explicitly, the pressure is given by 



(9) 



where 




P^ = -ip + pe-AB) (10) 



3 

is the pressure contribution of the quark matter, and 

P„ = (r„ - l)pe (11) 
where r„ is not necessarily equal to Fq, is that of the nuclear matter, and 

C (P - Ptr)/{pq - Ptr) iOT Ptr < P < Pq 

a = I (12) 
I 1 for Pq < p, 

is defined to be the scale factor of the mixed phase (Lin et al. 2006). We should notice that 
Pq is not in the usual form of MIT bag Pq = ^(ptot — 4P), where ptot is the (rest frame) 
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total energy density, and B is the bag constant. It is because in Newtonian simulation, 
we use the rest mass density p and specific internal energy e as fundamental variables in 
the hydrodynamics equations. The total energy density ptoh which includes the rest mass 
contribution, is decomposed as ptot = P + P^- We choose r„ < Fq in our simulations to 
take into account the possibility that the nuclear matter may not be stable during the phase 
transition process, and hence some quark seeds could appear inside the nuclear matter, or 
the convection mentioned in the early of this section may mix some quark matter with 
the nuclear matter. In the presence of the quark seeds in the nuclear matter, the effective 
adiab atic index will be reduced. The possible values of B ^^^ range from 145 MeV to 190 



MeV flDeGrand et al.lll975l : ISatzlll982 



Steffens et al.lll995l ). For p > pg, the quarks will be 



deconfined from nucleons. The value of Pg is model depen dent; it could range from 4 to 8 p„ 



(jCheng fc Dailll998bl : iHaensell 1200.4 iBombaci et al.ll2004l ). where p, 
the nuclear density. 



2.8 X lO^'' g cm ^ is 



In the simulations, we set Fq = 2, F„ = 1.85, S^/^ = 160 MeV, and pg = 9pnuc- The 
phase-transition density ptr ~ 2.1p„„c is defined to be the point where Pg is zero initially. 
This value is approximately equal to the value estimated in section 2.2. The total time span 
of each run is ~5 ms, the time step is 3.7 x 10""^ ms. The grid spacing is set to be dx = 0.28 
km and the outer boundary of the computational domain is at 27.5 km, which is about 
two times the stellar radius of our models. During the numerical simulation, a low-density 
atmosphere is added outside the neutron star for numerical stability purpose. The density 
and temperature of the atmosphere are 3 x 10^ g/cm^ and 0.003 MeV, respectively. 

Finally we want to remark that although eq.(9) is not the exact situation of quark 
matter distribution inside the star, it allows us to simulate the phase induced collapse and 
study the energy transport processes. The exact quark matter distribution can certainly 
affect evolution of star in detail and neutrino emission from the stellar surface quantitatively. 
However since eq.(9) should still represent the qualitative quark matter distribution therefore 
the simulated features produced by eq.(9) should be still qualitatively correct. 



3.2. Neutrino luminosity 



For a new born neutron star, the internal temperat ure is so hot that neu t rinos will be 
trapped inside the star for at least a few seconds (cf. IShapiro fc TeukolskyI (119831 ) for a 
general review). However, neutrinos very near the surface of the star can still escape from 
the star because the optical depth near the stellar surface is low. In fact we can define a 
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radius [Ru) called neutrinosphere, where the optical depth of electron neutrinos is given by 

drKeff{r) = l (13) 



eff 



where the effective optica l depth, Teff is defined as inverse mean free paths and the effective 
opacity, Kg// is given by (jjankall200ll ) 



(«:e//)(r) = 1.202 xlO-Vio(r) 



V4MeVy cm' 



(14) 



where pio{r) is the density in units of 10^° g/cm^, and Ty(r) is the temperature at r r espec- 
tively. The electron neutrino lum inosity emitted from the neutrinosphere is given by (jjanka 



19951 : iBalantekin fc Yukselll2005f ) 



16 



TrR^acT^, . 



(15) 



In our simplified calculation, we also assume equal luminosities for neutrino and antineutrino, 
hence the combined luminosity for a single neutrino flavor is 



Ll, + Ly 



-TrR^acTy 



(16) 



Since the neutrinosphere is very near the stellar surface, where the nuclear matter is 
described by an ideal gas EOS P„ = (r„— l)pe (see previous section), we thus approximate the 
temperature by T = 2'mBe/3k, where k is the Boltzmann's constant and rriB = 1-67 x 10~^^ 
g is the baryon mass. With the simulated density and temperature profiles at a given time, 
we obtain the value of each R^, from to 3 ms with trial-and-error method. Since the 
neutrinosphere is a function of both temperature and density, which oscillates with a period 
~ 0.32 ms; therefore it also oscillates with the same period. The temperatures and densities 
as a function of time are shown in Fig. 1 and Fig. 2 for the neutron star with 1.55 and 1.75 
Mq respectively. 

From these two figures it appears that the temperature oscillation and the density 
oscillation are almost close to 180° out of phase. Since the entire system is oscillating, 
therefore every quantity should oscillate with the same period. Also the equations governing 
the evolution of temperature and density are not identical, it is not surprising that they 
have a phase difference. However, the question is why they are almost 180° out of phase? 
This phenomenon can be explained qualitatively as follows. When the oscillatory cycle 
begins, matter is falling in until the end of the first half cycle. However, once the matter 
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is squeezed, the temperature is rising everywhere inside the star and the gradient of the 
temperature starts to drive the thermal energy outward. In fact most of thermal energy is 
generated in the core. Furthermore the definition of the neutrinosphere is that the product 
of the temperature and density is a constant. Therefore, when the density is decreasing with 
time, but the temperature is increasing with time in the first half of the cycle, there will 
always have a situation that the density is minimum and the temperature is maximum near 
the end of the first half of the cycle at the neutrinosphere. 

Using R„ and T,, at the neutrinosphere obtained from the numerical simulation, we 
compute the time evolution of the neutrino luminosity. The results arc shown in Fig. 3 and 
Fig. 4 for two different neutron star mass respectively. We can see that extremely intense 
neutrino pulses can occur when the temperature at the neutrinosphere is maximum while the 
density at the neutrinosphere is minimum. The pulsation period of the neutrino luminosity 
is equal to that of the temperature and density. In calculating the neutrino luminosity we 
have used the constraint L^, < Em where Em represents the rate of energy flow carried 
by fluid into the cell containing the neutrinosphere R„. Explicitly Em — ^TrR^ARr/St, 
where 6t is the dynamical time scale for the energy flow from one grid to another grid and 
T = O.bpv^ + pe is the total energy density at R^. 

The maximum neutrino luminosities range from 10^^ to 10^^ erg/s. The total energies 
carried away by neutrinos in the period of first 3ms are 6.75 x 10^^ ergs and 7.12 x 10^*^ ergs 
for the star with 1.55Mq and for the star with 1.75M0 respectively and the typical neutrino 
energy is > 30MeV during the peaks. 

3.3. Damping of the oscillation 

With the grid resolution (dx — 0.28 km) we used for the simulations, which is limited 
by the computational resource, we see that numerical damping becomes significant after 
about 3 ms. Recently Adikamalov et al (2008) have used a 2D general relativistic numer- 
ical code to study the gravitational-wave signals emitted from the phase-induced collapse 
of neutron stars and their results are very similar to our previous results (Lin et al. 2006). 
Since Adikamalov et al (2008) can achieve much higher resolution in 2D, they can perform 
longer timescale simulations with high accuracy. They conclude that the timescale of hydro- 
dynamical damping effects (e.g., due to mass-shedding) is typically a few tens to hundreds 
ms. 

However, there are still other physical damping effects which are not modelled in the 
simulations (Lin et al. 2006; Adikamalov et al 2008). Theoretically Wang & Lu (1984) first 
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pointed out that the dissipation due to nonleptonic reaction is of great importance and the 
stellar pulsations of the quark stars would be strongly damped via the following process 

s + u ^ u + d. (17) 



According to ISawyerl (119891 ) and iMadseru (119921 ). in the high-temperature limit, which 
is exactly our case, the bulk viscosity can be obtained analytically 



c 



(l - (1 - exp(-/3V2TV)) 



Uj2 + f]T^ ) 



a = 9.39 X lO^^m^/i^ g cm^'s 



-i„-i 



(19) 



/3 = 7.11xlO-V^(l + m^/4/i2)2 s-2 (20) 

where rris is the mass of strange quark mass in MeV, typically ranges from 100 to 300 MeV. 
fid is the down quark chemical potential, the typical value is around 235 MeV, assuming the 
nuclear matter density is 2.8 x 10^^ g cm~^. r is the perturbation period and u = 2tt/t. 



For the star with relative constant density, ISawyerl (119891 ) estimated the damping time 
of vibration: 

To = 30-Vi?\"^ (21) 

The average damping time is thus calculated with average temperature ((T)) and density 
((p)), which are given by. 

Jq^**^ Aiir'^T{r)p{r)dr 
Jg^"^ ATir"^ p{r)dr 
Jq^"^ ATTr'^p{r)dr 
J^'"- Anr^dr 

where Rtr is the radius where the matter transits to quark matter, i.e. p{Rtr) = Ptr- 

If we take ~ 140 MeV, (p) ~ 10"*^^ g/cm^ and (T) ~ 50 MeV, the damping time 
scale is ~ 10 s. This time scale is sufficient long enough to allow the neutrinos to carry away 
most of gravitational energy released from the phase-induced collapse neutron star {AEq ~ 
GM'^AR/R'^, where AR is the change of radius before and after the phase transition.). For 
the model presented in this paper, AR/R is ~ 0.2, which gives ~ 10^^ ergs, which does 
not conflict with the 1MB results. Furthermore, in our simulation neutrinos can carry away 
~ ]^g50 gj,gg 2 j-j^g assume that this is the only mechanism to damp out the oscillation 

energy, it only takes 3 ms x ^^"50 '^^^ ~ 3 s to take this amount of energy away by neutrinos 
with energy > 30 MeV. 



(T) = Z. : > (22) 



(p) = '° .K.. . . : > (23) 
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4. Cooling of quark stars and neutron stars 



The use of the thermally excited helical vortex waves that produce fast magnetosonic 
waves in the stellar crust, and which propagate toward the surface and transform into outgo- 
ing electromagne tic radiation, has allowed the direct determination of the core temperature 
of neutron stars (ISvidzinskyI |2003| ) . The core temperature of the Vela pulsar is T = 8 x 10^ 
K, while the core temperature of PSR B0656+14 and Geminga exceeds 2 x 10^ K respec- 
tively. These temperature estimates rule out the equations of state incorporating Bose 
condensations of pions or kaons and superfluid quark matter in these objects. Thermal X- 
ray radiation from neutron star soft X-ray transients in quiescence provides the strongest 
constraints on the cooling rates of neutron stars, and thus on the interio r composition and 
properties of matter in the cores of compact objects (IHeinke et al.ll2006l ). The analysis of 
the new (2006)and archival (2001) XMM-Newton observations of the accreting millisecond 
pulsar SAX J1808. 4-3658 in quiescence provides the most stringent constraints to date. Si- 
multaneous fitting of all available XMM data allows a constraint on the quiescent neutron 
star (0.01-10 keV) luminosity of L^s < 1-1 ^ 10'^^ erg/s. This limit excludes some current 



models of neutrino emission mediated by pion condensates (IHeinke et al.ll2006l ). and provides 
further evidence for additional cooling processes, such as neutrino emission via direct Urea 
processes, involving nucleons in the cores of massive neutron stars. Hence, these recent ob- 
servations show that the observed thermal lu minosity of some neutron stars does not ag ree 
with some proposed exotic cooling processes (ITsuruta Ill998l : lYakovlev et al.ll2003l . |2004| ). 



4.1. Energy dissipation mechanisms for strange stars 



Several physical processes that contribute to the energy emission from the bare quark 
star surface have been proposed. The most important of these energy loss mechanisms is 
the electron- positron pair creation, due to the strong electric field at the surfa ce of the 
stran ge star (Usoy 1998al lbl). Other s processes include neutral pion emission ( Ng et al. 
20031 ) . quark-quark bremsstr ahlung (IChmai et al.l Il99ll : ICheng fc Harkd l2003l ). the equi- 



l ibrium black body radiation ([chmai et al 



(jJaikumar et al.ll2004l : iHarko &: Cheng] |2005| ) . respectively. 



199ll ) and the electron-electron bremsstrahlung 



In the case of the energy loss via the production of e pairs, the energy fiux is given 
by F± = e±ri±, with e± = me + T and 



(24) 



1 /3 

where ep = (vr^ng) is the Fermi momentum of the electrons with number density rie, a is 
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the fine structure constant, ^ = 2^1^{eplT), J (0 = (1/3) In (1 + 2^-^) 1 (1 + 0.074£f + 

(vrf/G) (13.9 + ^)^, and Ar is the thickness of the electron layer f Usovll2001 ). iHarko fc Cheng 

(120061 ) have considered the mean value of the z-dependent chemical potential, Sp ^ (/ie (-z, T)) ^ 

1 

. For temperatures of around T = 10 MeV we 
3 MeV, a smaller value than the one considered in lUsovi (120011 ) . £i? ~ 18 



(l/d) ^/3/^{Vi/T) (l + 



obtain ep 

MeV. This choice of the Fermi momentum significantly reduces the energy flux from the 
electron-positron pair creation. 

Pions are created at the stellar surface due to the collision between the quarks and the 
bag, representin g the conversion of the quark kinetic energy to the mass of the pion cloud 
(INg et al.ll2003l ). The pions are assumed to decay via the processes 7r° — > 27 < — > e'^ + e" and 



TT 



6^ + 1^6 + 2z/a- The energy flux due to this process is given by 



where is the energy density of the pion field at the star's surface, which is fixed by the 
axial current conservation, and v-n-, the speed of the ern itted pions, is t;^ ~ A/2T/m^, where 
uIt, ~ 140 MeV is the mass of the pion (INg et al.ll2003l ). 



Electron-electron bremsstrahlung is an imp ortant energy loss mechanism for strange 
stars (jJaikumar et al.ll2004j : iHarko fc Cheng]|2005l ). For surface temperatures T < 10^ K, the 
photon flux exceeds that of the electron-positron pairs that are produced via the Schwinger 
mechanism in the presence of a strong electric field that binds electrons to the surface of 
the quark star. The average energy of photons emitted from the bremsstrahlung process can 
be 0.5 MeV or more, which is larger than that in electron-positron pair annihilation. The 
effect of the multiple and uncorrelated scattering on the radiation spectrum (the Landau- 
Pomeranchuk-Migdal effect ), together with the eff ect of the strong electric field at the surface 
of the star, was discussed in lHarko fc ChengI (120051 ). The presence of the electric field strongly 
influences the radiation spectrum emitted by the electrosphere. The radiation properties of 
the electrons in the electrosphere essentially depend on the value of the electric potential 
at the quark star surface. The effect of the multiple scattering, which strongly suppresses 
radiation emission, is important only for the dense layer of the electrosphere situated near 
the star's surface, and only for high values of the surface electric potential of the star. The 
bremsstrahlung emissivity of the electrosphere can be obtained as 



See) 
'Br 



75 
4^ 



ar^g Liy/a (r/) 



T 77 
In \ 

2m, 30 



7 Li7/2 iv) + 



rfLi„ (r/) 



dn 



n=7/2 



T\ (25) 



where g = 2is the statistical weight for electrons, rj = — exp [/Xg (T) /T], Li„ (r/) = J2T=i V^/^" 
is the polylogarithm function and 7 = 0.577216 is Euler's constant. 

The bremsstrahlung energy flux from the exterior electron layer of the strange star. 



coming out from a thin surface layer of thickness dz, is 



(ee) 



Br 



e^^l'dz/ir. Taking into 
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account the contribution of all layers we find 

fV {T) = - 4f T) dz = at) {T)T\ 



(26) 



where 



,(ee) 



(T) 



75 

47]-3 



Li7/2 [^(^,^)] 



In + ^-1] L17/2 (2, T)\ + 



2mp 



30 



(in 



n=7/2 



(27) 



The energy fiux from the quark-quark bremsstrahlung can be represented as -Fg-g = 
o'Br {nb,T)T'^, where aBr{nb,T) ^ g^Li"^ — exp (^{n'^nb)^^^ /T^ I {nb,T) X/ (27r)'^, A is the 
mean photon path in the quark matter, Uf, is the baryon number density at the surface of 
the star and 



a T 



1 + 4 (1 + 37r/4a,)' 



+ 



1 

r 



arctan2 (1 + Sn/Aas) — arctan 2x071^'^ 



1 + 37r/4a, 



an 



1/3 



+ 27 



D {I + M / Aas) - D (ran 



'I" 



(28) 



where is the strong coupling constant, a = In (l84Z /a^, r is the 

mean collision time between quarks, given by 



in 



8/3 



(29) 



which can be roughly estimated as r = l/nanv, an d the function D{x) is defined as 
D{x) = i [Li2 {—2ix) — Li2 {2ix)] (jcheng fc Harko 2003). ctq is the quark-quark elastic scat- 
tering cross section. Due to interference between amplitudes of nearby interactions, the 
bremsstrahlung emissivity from strange star surface is suppressed for frequencies smalle r 
than a critical frequency (the Landau-Pomeranchuk-Migdal effect) ( ICheng fc Harkd l2003l ) . 
The range of the suppressed frequencies is a function of the quark matter density at the star's 
surface and of the QCD coupling constant. For temperatures much smaller than the Fermi 
energy of the quarks the bremsstrahlung spectrum has the same temperature dependence as 
the equilibrium black body radiation. Multiple collisions could reduce the intensity of the 
bremsstrahlung radiation by an order of magnitude. The absorption in the semi-degenerate 
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electron gas can also significantly reduce the intensity of the quark-quark bremsstrahlung 
radiation and, consequently, the surface emissivity. The combined effects of multiple col- 
lisions and absorption in the electron layer could make the soft photon surface radiation 
of quark stars six ord ers of magnitude smaller than the equilibrium black body radiation 
(ICheng fc H^boosh . 



At temperature T strange matter is filled with electromagnetic waves in thermodynamic 
equilibrium with quarks. The quanta of the electromagnetic waves in plasma (transverse 
plasmons) have a characteristic dispersion relation oj j k) = + fc ^, where k is the wave- 

1 /2 

number and the plasma frequency tUp = e (Svmb/S/i) (iChmaj et al.lll991l ). Th e characteris- 



tic transverse plasmo n cutoff frequency can be estimated as ~ 20 — 25 MeV (jChmaj et al. 
199ll : iNg et al.ll2003l ). A beam of transverse plasmons hitting the edge of strange star mat- 
ter from inside will be partially reflected and partially will pass to the outer vacuum, being 
refracted. The energy flux of thermal equilibrium photons radiated from the bare strange 
surface is given by 



eq 



where g (cj) = (l/27r2) /J^^^ [l - [R^ + R{) /2] sin Q cos QdQ, with = sin^ {9 - 9o) / sin^ {9 + 9o 



Lu{LU^-Lul)giLu) [e-/^-l]"'ciw. 



(30) 



and R\\ = tan^ 



^^o) / tan^ {9 + 9o). 9q is deflned as 9o = arcsin 



spectively (jChmaj et al.lll99ll ). 



sin 9 



re- 



The variation of the energy fluxes F for these energy mechanisms are represented, as a 
function of the temperature, in Fig. 5. 

When the temperature of the quark star core drops below 10^ K, the strange mat- 
ter becomes superfluid. At this temperature quarks can form colored Cooper pairs near 
the Fermi surface and become superconducting. From the BCS theory it follows that the 
critical temperature at which the transition to the superconducting state takes place is 
Tc = A/1.76, where A is the pai ring gap energy (IBlaschke et al.l l2000l ) . An early estima- 
tion of A gave A ~ 0.1 — 1 MeV f Bailin &: Love 19841 ) . but some recent studies considering 
instanton-induced interacti ons between quarks estimated A ~ 100 MeV ( lAlford et al.lll999 



Rajagopal fc Wilczekll200lh . When the temperature of the star is below Tc, the emissivity 
is of the quark matter is modified by the superconducting effects. Since the collisions be- 
tween the qu arks and the ba g are suppressed, pion emissivity is suppressed by a factor of 
exp (— A/T) (jNg et al.ll2003l ). The same factor can be used to describe the suppression of 
the quark-quark bremsstrahlung due to the superconductivity of the strange matter. 

Fig. 6 shows the variation with temperature of the energy fluxes for various dissipation 
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mechanisms if quarks become superfluid. 



For T << ujp the equilibrium photon emissivity of strange matter is neghgible small, 
compared to the black body spectrum. Even for a small value of the gap pairing energy A 
the superfluidity of the quark matter strongly suppresses the quark-quark bremsstrahlung 
and pion emissivity of the star. 

Strange quark matter in the color-flavor locked (CFL) phase of QCD, which occurs for 
large gaps (A ~ 100 MeV), could be rigorously electrically neutral, d espite the unequal 



quark masses, even in the presence of the electron chemical potential (lAlford et al.l Il999 
Raiagopal fc Wilczekl [200ll ) . 



However, iPage fc Usovl (120021 ) pointed out that for sufficiently large the low density 
regime is rather expected to be in the 2-color-flavor Superconductor" phase in which only u 
and d quarks of two color are paired in single condensate while the ones of the third color 
and s quarks of all three colors are unpaired. In this phase, electrons are present. In other 
words, electrons may be absent in the core of strange stars but present, at least, near the 
surface where the density is lowest. Nevertheless, the presence of CFL effect can reduce 
the electron density at the surface and hence also signiflcantly reduces the bremsstrahlung 
emissivity of the electrons in the surface layer. 



4.2. Cooling of strange stars 

In this subsection, we consider the cooling of a quark star. For simplicity, we assume 
that the star is of uniform density and isothermal, that is, the core and surface temperatures 
are equal. The effect of the magnetic fleld on the cooling is also neglected. The thermal 
evolution of the quark star is determined by the equation 

dT " 

^^~dt ~ ~ ^ ~ ~ (■^PO'ir + Lqg -\- Lge + -Z^tt + -^eg + L^) , (31) 

i=l 

where Cy is the speciflc heat, T is the temperature, and Lpair is the electron-positron pair 
luminosity, Lqq is the quark-quark bremsstrahlung luminosity, Lee is the electron-electron 
bremsstrahlung luminosity, L^r is the pion emission luminosity, Leq is the luminosity due 
to the thermal equilibrium radiation, and Lj, is the neutrino luminosity, respectively. In 
the normal state of quark matter, the quark Fermi momentum ppg can be approximated 
as pfq = 235 (p/po)^^'^ MeV/c, where po is the nuclear density, and the speciflc heat o f 



the quark matter is given by Cq = 2.5 x 10^° (p/po)^^^ Tg erg cm ^ K ^ (llwamotol Il980l ). 



In the superfluid state, for 0.2Tc < T < Tc, the speciflc heat can be obtained as c, 
3.15c, (Te/T)exp(-1.76Te/T) [2.5 - 1.66 (T/Te) + 3.64 (T/Tj^j , while for T < 0.2Te the 
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specific heat is zero (see iNg et al.l (120031 ) and references tlierein). Tlie Goldstone excitations 
in tlie quark-gluon plasma contribute to tlie specific heat of the strange star, so that Cg_-y = 
3.0 X 10^^Ng_^Tg erg cm~^ K~^, where Ng_^ is the number of available massless gluon- 
photo n states, which are present even in the color superconducting phase (IBlaschke et al. 



20001 ). When the temperature is low, the specific heat of quark matter vanishes. However, 
the electrons are not affected, and at low temperatures their effect is import ant. The specifi c 



heat of the electrons is given by = 1.7 x lO^o {Yep/pof^ Tg erg cm K'^ flNg et al.ll2003l ). 



where ~ 



10 ^ is the electron fraction. 



Neutrinos are emitted by quark matter through the URCA process, d —>■ u + e + Ve 
u + e~ — * d + z/g-. The neutrino em issivity for thi s process can be obtained as Ed 



8.8 X 10^6 (p/po) Ye^^Tl erg cm ^ s ^ (jlwamotdll980l ) . In the superfiuid state the neutrino 



emissivity is suppressed by a factor of exp (— A/T). The exact form for the neutrino emis- 
sivity can be obtained as Ed = ('45 77r/840/i^°c^) Gjp cos^ 9c (1 — cosOue/a) PpuPpciPFe {ksT)^ 



erg cm ^ s ^ (IDuncan et al.lll983l ). where Gp is the Fermi weak coupling constant, 6c is the 



Cabibbo angle, a = (1 — 2ac/TT)~^^^ , Ppu, Ppd and Ppe are the Fermi momenta of the m, 
d quarks, and of the electrons, respectively. The angle between the u and e momenta, due 
is given by cos6'„e = ( P^ri ~ P^,, ~ -P p f^PpuPpe- The condition for the URCA process to 
occur is |cos6'„e| < 1 (INg et al.ll2003l ). In the case of the similar reactions for the s quarks, 
s — > u + e~ + i>e-, M + e~ — *• s + i^e-, the condition |cos6'^g| < 1, where the angle 6'^^ is defined 
as cos6'^g = {Pps — Ppu — Ppe) f^PpuPpe, is Satisfied only if the density of the quark matter 
is as high as 4 x 10^^ g cm~^. Therefore these reactions do not exist in the quark stars. 
The neutrino emissivity due to the photon gluon mixing, which can generated a massive 
photon-g luon excitation, is neg lected, since it is important only at temperatures higher than 
70 MeV (IBlaschke et al.l 120001 ). The cooling curve and the luminosity curve of the l.SSM© 
strange star are presented as the solid curve in Fig. 7 and Fig. 8 respectively. We can see 
the radiation luminosity at t ~ 20 years is below but with a high temperature 

T ~ lO^K and the spectrum is a thermal bremsstrahlung. 



4.3. Neutron star cooling process 



In this subsection we describe a simple cooling model of neutron stars including possible 
exotic matter core, e.g. meson condensate or quark core (see e.g. Cheng et al. 1992, 
Chong & Cheng 1993). We assume that the neutron star with l.SSM© and 10^ cm in 
radius is uniform density and has a isothermal core temperature (Tc). We adopt the simple 
relation between the cor e temperature and the surf ace temperature [Tg) given by the relation 
Tcs = 1.3{TgQ/ Qsu)^'*^^ (iGudmundsson et al.lll98ll ). The core temperature evolution of the 
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star is given by 



CydTc/ dt — —Ly — Lhiy{Ts 



The heat capacity Cy is given by (jMaxwelllll979l ) 
where 

C('=) = 1.9 X lO^^Mep^f Tg erg K"^ 
is the electron heat capacity with pi4 = g cm~^ and Tg = K, and 



(32) 

(33) 
(34) 



CniT) = 2.3 X IO^^Mq ^ ^ 
3.15C„(T) (f ) e-i-76T„/T 



Pil'^T^ erg K" 



2.5 



1.66^ + 3.64 , ^ 



for T>Tn 
erg for T <Tn 



(35) 

is the neutron heat capacity, with T„ the transition temperature of normal-superfluid neu- 
trons and m* the effective ma ss of neutrons. We adopted m*/m„ ~ 1 and T„ = 3.2 x 10^ K 
(ITakatsuka fc Tamagakilll97ll ). 

The major neutrino emissions are due to modified Urea process and NN-bremsstrahlung, 
where N is a ncleon, n or p. These mechanisms are relatively weak and the luminosity can 
be approximated by (Yakovlev & Haensel 2003) 



(36) 



where Lg is at the range of 4 x 10'^^ to 10*" erg s 



^40 



When considering some exotic neu trino processes, the em ission is greatly enhanced. The 
neutrino luminosity can be written as (lYakovlev et al.l 120031 ) 



(37) 



where Lf lies at the range 4 x 10*^ to 4 x 10*"^ erg s ^ for pion condensation process and 
4 X 10*^ to 4 X 10'^^ erg s~^ for kaon condensation and quark cooling process. 

The cooling curves and the radiation luminosities of neutron star with various cooling 
mechanisms are shown in Fig. 7 and Fig. 8 respectively. We can see that the neutron star 
without exotic cooling processes, e.g. pion condensation or kaon condensation, their thermal 
luminosity is much larger than 10^* erg/s. But the the exotic cooling processes can cool the 
star much faster and the thermal luminosity can be much lower than 10^*erg/s. On the other 
hand, the surface temperature of the neutron star with exotic cooling processes is very low 
< lO^K and it is thermal spectrum. 
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Discussions and final remarks 



The detection of neutrinos from SN1987A lias suggested that a compact object should 
be formed at the center of SN1987A. If the time delay and energy distribution difference 
between neutrinos detected by KAM-II and 1MB are genuine together with the absence of 
hot compact object at the center of 1 98 7A, which implies a rapidly cooling cold compact, we 
have proposed that the neutrinos detected by 1MB may result from a delay phase-induced 
collapse after the neutron star is formed. The final compact object may be a neutron 
star with exotic matter core or a strange star, either of them is a rapidly cooling object. 
Consequently the current radiation luminosity from SN1987A should be much lower than 
Iq3a QYg/s. However, we have argued that the phase transition from neutron star to strange 
star may be the most possible phase transition occurred in 1987A. Such collapse process 
can generate temperature and density oscillation on the neutrinosphere. By using a three 
dimensional Newtonian code with a toy model, in which the phase transition is taken to 
be instantaneous and the equation of state is assumed to be a mixed phase, we find that 
intense pulsating neutrinos with energy >30MeV can be emitted. This phenomenon should 
be generic because the compression generates more heat in the core and the heat flow is 
outgoing during the contraction phase, whereas the matter is flowing in. That should create 
a situation that when the temperature is maximum and the density is minimum at the 
neutrinosphere. It is very interesting to detect intense pulsating hot neutrino emission from 
supernova explosions occurrin g from nearby galaxies. The future neutrinos experiments, like 
ANTARES, and IceCube (see iHalzenl (120061 ) for a review) may be able to detect such signals. 



In calculating the cooling of strange star, we have assumed a bare strange star, i.e. a 
strange star without crust. This assumption may be justified by the following reasons: On 
the top of the quark matter, a thin cr ust can exist as long as the electron density within it is 
smaller than that in the quark matter (lAlcock et al.lll986l ). Such a baryonic crust is, however, 
much thinner than the neutron star crust. The mass of the crust on the quark star is at 
most 1O~^M0, because the density of the crust bottom cannot exceed neutron drip density, 
and it may extend to up to 250 — 300 m. The normal crust may occur, for instance, due to 
the accretion of normal matter onto a bare strange star. The oscillation spectr um of strange 
stars with crust differs from the spectrum of neutron stars (jChuganovl l2006l ) . If detected, 
acoustic oscillations would allow one to discriminate between strange stars with crust and 
neutron stars, and constrain the mass and radius of the star. The recent detection of seismic 
vibrations in the aftermath of giant flares from two magnetars (highly magnetized compact 
stars) is a major breakthrough in observational astronomy. The oscillations excited seem 
likely to involve the stellar crust, the properties of which differ dramatically for strange stars. 
The resulting mode frequencies for st range stars cannot be reconciled with the observations, 
for reasonable magnetar parameters (IWatts fc Reddyl 120071 ). On the other hand, young and 
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hot strange stars are the source of a powerful pair wind, consisting of electron-positron pair s 



and photons, created by the Coulomb barrier at the quark surface (lAksenov et al.ll2005l ). 
Photons dominate in the emerging emission, and the emerging photon spectrum is rather 
hard and differs substantially from the thermal spectrum expected from a neutron star with 
the same luminosity. The total luminosity in the case of a stationary, spherically outflowing, 
pair winds is in the range of L = 10^^ — 10^^ ergs/s. These results have direct relevance 
to the emission from hot, bare, strange stars. For L > 2 x 10^^ erg/s, photons dominate 
the emerging emission. As L increases from 10^^ to 10^^ ergs/s, the mean photon energy 
decreases from ~ 400 — 500 keV to 40 keV, while the spectrum changes in shape from a 
wide annihilation line to being nearly blackbody with a high energy (> 100 keV) tail. Such 
a correlation of the photon spectrum with the luminosity, together with the fact that super- 
Eddington luminosi ties can be achieved , might be a good observational signature of hot. 



bare, strange stars (lAksenov et al.ll2003l ). At the moment of its formation, a strange star is 
very hot. The temperature of the interior may be as high as a few of 10^^ K. The neutrino 
luminosity of the young quark star is of the order of 10^^ erg/s, the rate of mass ejection from 
such a hot compact object is very high, and the normal matter envelope is blown away by 
radiation pressure in a few seconds. High temperatures also lead to a considerable reduction 
of the Coulomb barrier, increasing the tunneling of nuclei through the barrier toward the 
surface. Therefore, it is nat ural to exp ect that the surface of a young strange star to be 



nearly (or completely) bare (IUsovll200ll ). Since the strange quark matter at the surface of 
the star is bound via strong interactions rather than gravity, such a star can radiate at the 
at the luminosity greatly exceeding the Eddington limit. Therefore, due to the radiation 
pressure, a normal matter crust cannot be built around a young strange star. 

We have also assumed that the phase transition should take place after the neutron 
star is formed. We argue that the delay time scale of the phase transition is of order of 
the neutrino trapping time scale based on the following reasons: The theory of type II 
supernovae predicts the gravitational collapse of a degenerate core on timescales of the order 
of tcoii ~ (Gp)"'/' ~ 4 X 10- '^Pi2-5 ~ 0.1 s, followed by the formation of a hot [ksT ~ 50 
MeV) protoneutron star, with a radius of around 50 km, composed o f Nt ~ 10^^ baryons. 



and with the electron lepton number Lg ~ 0.35iVb (iHaensel et al.l 120071 ). The initial mass of 
the collapsing core is Mcore ~ A'tmo, where ttlq = 1.6586 x 10~^^ g is the mass of the ^^Fe 
nucleus, divided by 56. The gravitational mass of the protoneutron star is only slightly lower 
than Mcore, because the gain in the gravitational energy is compensated by the increase of 
the internal energy contained mostly in the strongly degenerate neutrinos z/g, trapped in the 
stellar interior. The diffusion time-scale of the neutrinos may be estimated by assuming that 
coherent scattering is the dominant opacity source, so that t^iff ~ ^A^Ngcatt/c, where A™'' 
is the mean free path in the sea of heavy nuclei {A, Z) and Ngcatt > 1 is the number of 



-27- 



scatterings experienced by the neutrino prior to escape. By taking into account the exphcit 
expressions of the neutrino mean free path and of the number of scatterings, the neutrino 
diffusion time can be written as tdiff ~ 0.08pi2 s. At sufficiently high densities tdifj » tcou 
(IShapiro &: TeukolskyI Il983l ). The two timescales becomes comparable at t^iff ~ tcou, a 



condition which is satisfied at a density ptrap ~ 1-4 x 10^^ g/cm^. 

Neutrino trapping has enormous implications for the core collapse. For p > ptrap, most 
of the neutrinos from electron capture remain in the matter, and the lepton number per 
baryon does not change. The neutrino distributions approaches an equilibrium Fermi-Dirac 
distribution. The neutrino luminosities are greatly reduced by trapping. Once the center of 
the core exceeds nuclear densities of the order of pnud = 2.8 x 10^"^ g/cm^, thermal pressure 
and nuclear forces cause the equation of state of the matter to stiffen, preventing further 
collapse. Most of the gravitational binding energy of the core is released in the form of 
neutrinos, following the collapse to nuclear densities. In the absence of neutrino trapping, 
the total binding energy would be completely emitted as neutrinos in a collapse timescale, 
the time for the core to contract from 2Rnud to Rnud, where Rnud ~ 12 km for M ~ 
Mq. Accordingly, the neutrino luminosit y would then achieve i ts ma ximum possible value. 



Lu,ma.x = GM"^ / Rnud/ 1 coll ~ 10^'' crg/s (IShapiro fc TeukolskyI Il983l ). In reality, neutrino 



trapping forces the liberated gravitational potential energy to be emitted on a much longer 
diffusion timescale, tdiff » tcoii at p ~ Pnud- As a result, the actual neutrino luminosity is 
of the order L^, = GM"^ /Rnud /tdiff ~ 10^^ erg/s. Thus during the advanced stages of the 
collapse, the neutrinos are unable to stream freely out of the core. The bulk of the liberated 
gravitational energy must therefore be converted into other forms of internal energy (e.g. 
thermal energy, energy of the excited nuclear states, bounce kinetic energy etc.), rather than 
being released immediately in the form of escaping neutrinos. 

In order for the phase transition to take place in the neutron star, the density inside 
the core must satisfy the condition p > ptrans , as discussed in Section 2. Since the core 
temperature during the neutrino trapping period is very high, the thermal pressure will 
give a significant contribution to the total pressure, thus lowering the central density of the 
neutron star, so that p < ptrans- In order to satisfy the thermodynamical conditions for the 
phase transition, the core must cool down so that the core density increases. To estimate the 
density increase resulting from the temperature decrease in the core of the neutron star we 
solve the TOV equation by using a temperature dependent EOS, i.e. Shen et al. (1998), to 
obtain core density as function of temperature for a fixed stellar mass. In implementing this 
EOS, we have added an ideal electron gas contribution. We have also imposed the condition 
of beta equilibiurm to determine the electron fraction. The left panel of Fig. 9 shows the 
energy density profiles for a star model with M = l.SSM© at different temperatures. Note 
that the temperature profiles are determined by the general relativistic isothermal condition. 
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The right panel of Fig.9 shows the results for a neutron star with M = 1.7Mq. In Fig. 9, we 
can see that indeed when core temperature decreases by a factor of 2, the core energy density 
can increase by more than 20%. We can see that these changes are sufficiently large to trigger 
the phase transition when the core temperature decreases after the neutrino trapping time 
scale. 

We conclude that the abnormally high energy neutrinos and the time delay detected by 
the 1MB detector can be explained by the neutrino trapping in the newly formed very hot 
neutron star, followed by a phase transition. After phase transition if the compact object is 
a neutron star with exotic matter core, its surface luminosity is lower than lO^^erg/s and it 
should have black body thermal spectrum with temperature < lO^K. On the other hand if 
it is a bare strange star, its surface luminosity is also lower than lO^^erg/s but the radiation 
spectrum is a thermal bremsstrahlung spectrum with temperature > lO^K. However, the 
nature of the phase transition cannot be inferred from the present data, since the observed 
features cannot differentiate or discriminate between the different phase transition models. 
We suggest that the future observations on spectral features of the compact remnant of 
1987A could provide better clues to elucidate the nature of the phase transition that might 
take place in SN1987A. 
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Fig. 1. — Left: Neutrinosphere temperature versus time. Right: Neutrinosphere density 
versus time. The initial mass of neutron is 1.55M0 
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Fig. 2. — Left: Neutrinosphere temperature versus time. Right: Neutrinosphere density 
versus time. The initial mass of neutron is l.TM© 
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Fig. 4. — Neutrino luminosity as a function of time of a star with l.TM©. 
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Fig. 5. — Temperature variation of the energy fluxes (in a logarithmic scale), emitted by the 
bare quark star surface, via different radiation mechanisms: electron-electron bremsstrahlung 
flux F^^^ (solid curve), electron-positron pair creation energy flux F± (dotted curve), quark- 
quark bremsstrahlung energy flux Fq^g (dashed curve), energy flux F.^ due to pion emission 
(long dashed curve) and black body radiation energy flux F^f, (ultra-long dashed curve). 
For the surface electrostatic potential of the quark star wc have chosen the typical value 
Vj = 14 MeV, while the thickness of the electron layer was taken d = 1000 fm. For the 
electron-positron energy flux the Fermi energy of the electrons Sp — 18 MeV=constant. For 
the energy density of the pion field at the strange star surface we have adopted the value 
^ 3.417 X 10^ MeV^ ^ 7.1 x 10=^^ erg/cm^. 
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Fig. 6. — Temperature variation of the energy fluxes (in a logarithmic scale), emitted by 
the superfluid bare quark star surface, via different radiation mechanisms: electron-electron 
bremsstrahlung flux F^^^^^ (solid curve), electron-positron pair creation energy flux {F±) 
(dotted curve), quark-quark bremsstrahlung energy flux F^^^^^ (dashed curve), energy flux 
F^^'^P'' due to pion emission (long dashed curve) and the thermal photon equilibrium radiation 
energy flux Fgg (ultra-long dashed curve). For the surface electrostatic potential of the quark 
star we have chosen the typical value V thickness of the electron layer was taken d — 1000 
fm. The energy gap A = 1 MeV. 
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Fig. 7. — Surface temperature as a function of time. The heavy sohd hne is the bare strange 
star coohng curve. The hght sohd Unc, the dotted hne and the dashed hne are the coohng 
curves of neutron star with modified URCA, kaon condensate and pion condensate respective. 
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Fig. 8. — The surface photon luminosity as a function of time. The heavy sohd hne is 
the bare strange star coohng curve. The hght sohd hne, the dotted hne and the dashed 
hne are the coohng curves of neutron star with modified URCA, kaon condensate and pion 
condensate respective. 
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Fig. 9. — Left: The energy density profile of stars with mass 1.55 Mq at different surface 
temperatures Ts = 1, 10, 15, 19 MeV. The central temperatures of the stars are Tc — 
1.49, 14.7, 21.5, 26.0 MeV respectively. Right: The energy density profile of stars with mass 
1.7 Mq at different surface temperatures Tg = 1, 10, 15, 20 MeV. The central temperatures 
of the stars are = 1.57, 15.5, 22.8, 28.9MeV respectively. 



